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We here argue that the “knee” of the cosmic ray energy distribution at E c rts 1 PeV represents 
a second order phase transition of cosmic proportions. The discontinuity of the heat capacity per 
cosmic ray particle is given by Ac = 0.450196 ks- However the idea of a deeper critical point 
singularity cannot be ruled out by present accuracy in neither theory nor experiment. The quantum 
phase transition consists of cosmic rays dominated by bosons for the low temperature phase E < E c 
and dominated by fermions for high temperature phase E > E c . The low temperature phase arises 
from those nuclei described by the usual and conventional collective boson models of nuclear physics. 
The high temperature phase is dominated by protons. The transition energy E c may be estimated 
in terms of the photo-disintegration of nuclei. 

PACS numbers: 13.85.Tp, 96.50.S, 64.60.Bd 


I. INTRODUCTION 

We have recently [lj discussed the power law energy 
exponent values {a} in the cosmic ray particle flux 
distribution [I, [3j 
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The more detailed measured exponent a in reality takes 
on two values depending upon the cosmic ray particle 
energy. For energies lower than the crossing energy E c ss 
1 PeV, the exponent takes on a boson value of 2.701178. 
For energies higher than the crossing energy the exponent 
takes on the fermion value of 3.151374. The crossing 
energy is, of course, the location of the so called “knee” 
in the energy distribution. 

In that the energy distribution depends on the heat 
and thereby entropy of evaporation of bosons from the 
cosmic ray source, there exists a quantum phase transi¬ 
tion of cosmic proportions at the crossing energy. It is 
evidently a quantum phase transition since the order pa¬ 
rameter involves the difference between Bose and Fermi 
statistical phases. In virtue of the experimental continu¬ 
ity of the energy distribution and thereby the entropy, 
the phase transition is higher than first order. To the 
theoretical and experimental present accuracy, the phase 
transition is of second order with discontinuities in the 
second derivatives of the entropy although more compli¬ 
cated non-analytic singularities cannot be ruled out. In 
Sec HU the thermodynamic properties of ultra-relativistic 
ideal gases are reviewed. The heat capacity discontinuity 
is described in Sec lII Al 

The boson phase arises from those evaporation nuclei 
described by the conventional collective Boson models of 
nuclear physics (4j. There exist pairing correlations in odd 
odd nuclei made up of deuterons. Correlations between 
two spin one deuterons lead to spin zero alpha parti¬ 
cles and so forth all within the pairing condensate. The 
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condensate resides near the surface of evaporating high 
baryon number A> 1 nuclei. For example, a neutron 
star itself is merely a nucleus of extremely high baryon 
number (A 1) with superfluidity (and supercon¬ 

ductivity) in the neighborhood of the nuclear surface^. 

The “partons” from the pairing condensate are evi¬ 
dently the nucleons. The partons then turn into fermions 
for the phase that exists above the crossing energy. For 
energies above the knee, the cosmic rays must be com¬ 
posed mainly of protons. In order to comprehend the 
phase transition, it is necessary to understand how pho¬ 
tons can photo-disintegrate the compound boson odd odd 
nuclei ultimately into its nucleon parts[c*]. In Sec lIII Al 
we employ a simple physical kinetics model of photo¬ 
disintegration to estimate the crossing energy. There is 
satisfactory agreement with experiment. A general view 
of the quantum cosmic ray phase transition is given in 
the concluding SeclIVI 


II. IDEAL GAS THERMODYNAMICS 

In the ultra-relativistic limit wherein the single particle 
energies are much larger than either the rest mass energy 
and the source environmental chemical potentials, the 
energy £, pressure P and volume V of the gas are related 
by 


£ = 3 PV. (2) 

Eq.© holds strictly true for “massless” gas at “zero” 
chemical potential. For such a case, the equation of state 
for the ideal gas reads 


PV = Nks'd wherein d = T , (3) 


wherein 

cxmb = 3 (Maxwell — Boltzmann Statistics), 
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asE ~ 2.701178 (Bose — Einstein Statistics), 
cvfd ~ 3.151374 (Fermi — Dirac Statistics), (4) 


with the solution 


i.e. 


(Maxwell — Boltzmann Statistics) => 

PV = J\fk B T, 

(Bose — Einstein Statistics) =>- 

PV « 0.9003926 J\fk B T, 
(Fermi — Dirac Statistics) =>• 

PV« 1.050458 Nk B T, 


( 5 ) 


wherein Af is the number of particles. One may verify 
the Bose-Einstein case by computing the pressure of a 
photon gas, i.e. black body radiation. One may verify 
the Fermi-Dirac case by computing the pressure of a gas 
of Weyl neutrinos. 

The pressure from the Bose-Einstein gas is slightly 
lower than that of a Maxwell Boltzmann gas in that the 
quantum statistics describes an attraction. The pressure 
from the Fermi-Dirac gas is slightly higher than that of 
a Maxwell Boltzmann gas in that the quantum statistics 
describes a Pauli exclusion repulsion. Although the value 
change in alpha due to quantum statistics is small, 


Aa = ( olfd — oi B e) ~ 0.450196, 


( 6 ) 


Eq. © is entirely responsible for the quantum cosmic ray 
second order phase transition. 


A. Heat Capacity 


The heat capacity at constant volume obeys 


Ca ^ = ( w 
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In virtue of Eqs.® and (HI in a regime in which a is 
uniform, the heat capacity per particle c = (Cuy/Af) is 
related to the energy per particle E = {£/AT) 

E = cT = ak B T. (8) 

The heat capacity discontinuity is thereby 

Ac = k B Aa « 0.450196 k B (9) 

in virue of Eq. © • 

B. Entropy 

If s(E) represents the energy per particle, then 
(1/T) = ( ds/dE ) determines the temperature. Eq.© 
then turns into a differential equation; It is 
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s(E) = ak B In ( — ) . 


( 11 ) 


In a single phase regime, the entropy of evaporation de¬ 
termines the energy distribution [ij via 
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On either side of the crossover energy, the distribution of 
energy is determined by Eq. (fl2l) . 


III. PHOTO-DISINTEGRATION 

The cross-over energy will here be estimated on the 
basis of the physical kinetics. From this point of 
view, one goes from the lower temperature boson phase 
to the higher temperature fermion phase via photo¬ 
disintegration processes. 

A. Relativistic Kinematics 

Consider a process wherein an initial compound nu¬ 
cleus I is hit by a photon 7 and disintegrates into final 
fragments F, 


7 + 1—if, 

with the four momentum conservation 

hkry + pi = PF- 


(13) 


(14) 


The invariant mass squared s of the reaction Ea. (fl3l) 
obeys 

c 2 s = —pp = —( pi + hk 7 ) 2 = — p 2 — 2 Tikry ■ pi (15) 
for a real (not virtual) photon k 2 =0. Thus 


c 2 (Mp — M 2 ) = —2 hk-y ■pi . 


(16) 


In the rest frame of the initial compound nucleus, the 
photon frequency w 1 obeys —k 7 -pi = Mjw-., so that 


hu7 7 = (Mp — Mi)c 2 


Alp + AIi 
2 Wi 


(17) 


wherein the term in square brackets on the right hand 
side of Eq. HU takes into account the total recoil of the 
final state fragments. 

In the rest frame of the cosmic ray source, the photon 
frequency is given by 


y (l — cos#) = 


c 2 Ali 

E 


ftw 7 (l — COS0) = c 4 
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E = ak B 


(18) 
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wherein E is the energy of the initial compound nucleus 
and 9 is the angle between the photon and compound 
nuclear three momenta. Averaging over angle and types 
of final fragmentation products yields 0 


htu^ = c 4 


Mj. - Mf 


2 E 


(19) 


The central Eq. (THU) of this section enables us to estimate 
the crossing energy. 


B. Crossing Energy 


We have discussed above that ~ 3fcsT 7 , wherein 
T 7 the radiation temperature of the source. The crossing 
energy can thereby be estimated as 


E c 


.4 Mj. - Mf 

6k 


( 20 ) 


The orders of magnitude involved in the estimate Eci. (l2(l 
are 


E r 


(100 MeV) 1 2 3 4 5 
10 eV 


1 PeV. 


( 21 ) 


The crossing energy in Ea. (l2T1) is in satisfactory agree¬ 
ment with experiment. 


transition from a lower temperature boson dominated 
cosmic ray beam to a higher temperature fermion dom¬ 
inated cosmic ray beam. The lower energy boson dom¬ 
inated regime owes its existence to the collective boson 
nuclei at the sources. These lead to the quantum sym¬ 
metry energy contribution to the semi-empirical nuclear 
mass formula [7]. The higher energy fermion phase owes 
its existence to the free single nucleons which are the par- 
tons from which the bosonic nuclei were constructed. As 
the energy per particle increases there is a de-confinement 
quantum statistical phase transition as the bosons are 
disintegrated into constituent fermions. 


There are many important applications of this picture 
of cosmic ray structure. Some light has been shed upon 
the structure of ordinary nuclei in terms of collective bo¬ 
son models. The notion of a superconducting shell near 
the surface of neutron stars has been explored. Clearly, 
nuclear transmutations and the cosmic ray sources of dif¬ 
ferent chemical elements should now be studied anew. 
With regard to these future developments, we are at the 
beginnings. However, we hope that our studies of cos¬ 
mic ray dynamics have clarified their deep connections 
to fundamental processes through rather precise deter¬ 
minations of the two correct critical indices along with 
an estimate of the cross-over energy . 
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